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Abstract 

Magnetic activity in M-dwarfs present enigmatic ques- 
tions: On the one hand they have higher field strengths 
and larger filling factors than the magnetic field on the 
Sun, on the other hand, they are fully convective and 
their atmospheres are more neutral, hence they do not 
have an undershoot layer for magnetic flux storage and 
as we show here, cannot have small-scale dynamo action 
in their photospheres either. We present a discussion of 
these facts and propose a new numerical model to inves- 
tigate M-dwarf magnetism. 
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1. Introduction 

M-type dwarf stars display a high degree of magnetic ac- 
tivity. Both indirect tracers as X-ray and H a emission, as 
well as Zeeman broadening of magnetically sensitive pho- 
tospheric lines has been observed (Johns-Krull & Valcnti 
1996, Kochukhov et al. 2001). The photospheric lines show 
no rotational modulation or net polarization, indicating 
that fields of small scale relative to the stellar radius ex- 
ist on the surfaces of M-dwarfs. These fields can have 
strengths of a few kG and can cover a substantial area 
with filling factors up to 50%. 

On the theoretical side of things, it has long been 
considered a problem that M-dwarfs are fully convective, 
lending no helping hand to the storage of magnetic flux 
through the presence of a stable stratified convective un- 
dershoot layer (e.g. van Ballegooiien 11982)) . Furthermore, 
as we show below (see also Dorch & Ludwig :2Q02|, it seems 
that their upper atmospheres are not capable of providing 
a small-scale network magnetic field, similar to that in the 
Sun, sometimes assumed to be stemming from local small- 
scale dynamo action in the photosphere (Cattaneo^999). 
However, recent global scale simulations of M-dwarfs seem 
to render the storage problem obsolete (Dobler 2004 ) , and 
in this contribution we propose future possibilities regard- 
ing numerical simulations of M-dwarf magnetic activity. 
First, however, we discuss results from models of M-dwarf 
local dynamo action. 
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Figure 1. Snapshot of the vertical component of the mag- 
netic field at optical depth unity (contours), emergent in- 
tensity (grey tones), and the velocity field (arrows) to- 
wards the end of an evolutionary seguence of a model with 
Rm=20. 



2. Local small-scale dynamo action 

In order to study local small-scale dynamo action in the 
photosphere of M-dwarfs we adopted a velocity field from 
a radiation-hydrodynamics simulation of a prototypical 
M-dwarf atmosphere from Ludwig et al. (2002), see Fig. 
^ Here T c s — 2800 K, log g — 5.0 and the chemical com- 
position is solar. Magnetic fields can potentially influence 
the flow structure significantly provided the field is strong. 
However, this corresponds to the full MHD case which is 
beyond the scope of this contribution. We selected a se- 
quence which comprises roughly 10 convective turn-over 
time scales. The sequence consists of 150 snapshots of the 
flow field, each comprising 125 x 125 x 82 grid points cor- 
responding to 250 x 250 x 87 km 3 . At any instant in time 
about 10 granular cells were present in the computational 
domain, ensuring a statistically representative ensemble. 

To compute the evolution of the magnetic field in re- 
sponds to the adopted velocity field, we assume the kine- 
matic regime of MHD, where one can neglect the back- 
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Figure 2. The magnetic Reynolds number as a function of 
Rosseland optical depth in a T c g = 2800K M-dwarf (solid 
line) and a solar model atmosphere (dashed line). Note the 
low magnetic Reynolds number in the M-dwarf model, pri- 
marily reflecting the low electric conductivity of the stellar 
gas in the rather cool M-dwarf atmosphere. 

reaction of the magnetic field on the fluid motions. Then 
solving the MHD equations reduces to the problem of seek- 
ing the solution to the time-dependent induction equation: 

- = Vx(uxB) + ,V 2 B, (1) 
at 

where u is the prescribed velocity field, B the magnetic 
field, and rj the magnetic diffusivity that. We assume a 
spatially constant magnetic diffusivity that we vary by set- 
ting the magnetic Reynolds number Rm defined as Rm = 
U£/rj (where U is a characteristic velocity, I a character- 
istic length scale, and r\ the magnetic diffusivity). 

We solve Eq. using staggered variables on the grid 
of the hydrodynamical flow field. The numerical method 
was originally developed by Galsgaard and others (Gals- 
gaard & Nordmnd ll997[l for general MHD purposes, but 
a special version is the code used by Archontis, Dorch, 
& Nordlund |2003| for studying dynamo action in pre- 
scribed flows. The kinematic approximation to the MHD 
equations becomes inaccurate when the magnetic field be- 
comes sufficiently strong, limiting us to weak fields. 

In the following we discuss kinematic MHD models 
with Rm as low as 20. In the astrophysical context, low 
values of Rm are uncommon, mostly due to the large spa- 
tial scales usually involved. However, in M-dwarf atmo- 
spheres we can be confronted with the situation of rather 
low Rm (Meyer & Meyer-Hofmeister 1999): Fig. [21 shows 
Rm in the M-dwarf model (£ = 80 km, v = 0.16 km/s) 
as well as the Sun (I — 1500 km, v = 2.4 km/s). Rm pri- 
marily reflects the run of the electric conductivity in the 
atmosphere, which in turn is mostly controlled by the elec- 
tron to gas pressure. The conductivity has been evaluated 
assuming a weakly ionized plasma. 
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FigureS. The total magnetic energy Em/Eo as a function 
of time in units of 150 seconds (the typical convective 
turn-over time). Five models are presented, with differ- 
ent amounts of diffusion: Rm = 800 (full), 400 (dotted), 
200 (dashed), 100 (dashed dotted), and 20 (dashed triple 
dotted curve). 

At sufficiently cool temperatures Rm reaches order 
unity in the surface layers. This is a consequence of the 
declining electron density, the shrinking of spatial scales, 
and smaller convective velocities: This refers to the sur- 
face layers, but qualitatively we expect a strong increase 
of Rm with depth, and beyond a certain depth the regime 
of Rm > 1 is reached again. However, gas motions in this 
depth will generally be slower and the tangling of magnetic 
field lines less rapid, which may reduce the efficiency of 
chromosphcric and coronal heating. Whether this plays a 
role for the observed decline of stellar activity at the tran- 
sition from M- to L-dwarfs (Gizis et al. 2000) is presently 
a matter of debate (Mohanty et al. 2002, Berger 2002). 

A flow is a fast kinematic dynamo when the exponen- 
tial growth rate 7 is positive. That is, fast dynamo action 
requires a continuous increase of magnetic energy, even 
in the limit of vanishing diffusivity. This limit is relevant 
because most astrophysical systems have Rm 3> 1 and 
small but non-zero 77. It is believed that turbulent astro- 
physical systems are fast dynamos operating at Rm 3> 1. 
When Rm increases, the length scale of magnetic islands 
decreases and scales as Rm~ 5 . There is a maximum value 
of Rm, which can be achieved with our numerical reso- 
lution Ax = 2 km, and therefore the largest magnetic 
Reynolds number that we can allow is of the order of 400 
corresponding to the Nyquist wavelength 2 Ax. In the fol- 
lowing, we concentrate on situations with Rm < 400. 

In the absence of non-linear effects, in case of a dy- 
namo, one expects a continuing exponential growth of 
magnetic energy Em- Figure [3] shows the result in terms of 
Em(£) for five different models corresponding to varying 
Rm: Most of the models in fact are dominated by decaying 
modes, with negative growth rates. In particular, the case 
with Rm = 20 is clearly an example of an anti-dynamo: 
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Figure 4- Illustration of loops formed in a solar corona 
simulation using a synthetic TRACE 195 A filter. 

The diffusion works faster than the flows can sweep up the 
field and concentrate it in the inter-granular down-draft 
lanes, and the dominant magnetic mode is a decaying one. 
fncreasing Rm decreases 7, so that for Rm = 200 the de- 
cay time is 40 time longer than at Rm = 20. fn terms 
of providing dynamo action, the most promising cases are 
those with lower diffusion and Rm > 300-400: At Rm = 
400 a growing mode seems to be dominating. 

In the high diffusion case the magnetic field varies 
smoothly across the domain, and its power peaks on scales 
approximately £ ss 50 km when considering B z (see Fig. 

At higher Rm more small scales are generated mostly 
around the down-draft lanes. 

3. Discussion and conclusion 

Little is known about the structure of magnetic fields 
in the photospheres of M-dwarfs. From an observational 
point of view one would like to get some input from the- 
ory that would alleviate the problem of disentangling field 
strength and filling factor. To this end we started to inves- 
tigate the kinematic effect of the convective velocity field 
on a magnetic seed field. This has become possible due to 
recent progress in the hydrodynamical modelling of atmo- 
spheres in the regime of cool M-dwarfs (Ludwig, Allard, 
& Hauschildt 2002). 

We performed kinematic studies of the evolution of 
small-scale magnetic fields in the surface layers of M-dwarfs 
We solved the induction equation for a prescribed veloc- 
ity field, magnetic Reynolds number Rm, and boundary 
conditions in a Cartesian box, representing a volume com- 
prising the optically thin stellar atmosphere and the up- 
permost part of the optically thick convective envelope. 
We find dynamo action for Rm > 400 and growth time 
scales of the magnetic field are comparable to the convec- 
tive turn-over time scale (~ 150sec). The convective ve- 
locity field concentrates the magnetic field in sheets and 
tubular structures in the inter-granular down- flows. Per- 
haps surprisingly, Rm is of order unity in the surface layers 
of cooler M-dwarfs, rendering the dynamo inoperative. In 



all studied cases we find a rather low spatial filling factor 
of the magnetic field. 

— The geometry of small-scale magnetic fields looks sim- 
ilar to the situation for the Sun. The basic reason is 
that M-dwarf granulation is qualitatively similar to so- 
lar granulation. 

— There are differences due to potentially quite different 
magnetic diffusivities. 

— Depending on Rm, we get or do not get local dynamo 
action. 

Future work along the lines of understanding M-dwarf 
magnetism will involve a more complicated treatment than 
that discussed above: It is our goal to combine both nu- 
merical theoretical modelling and detailed observations to 
begin to understand how M-dwarfs dynamos operate. We 
will use numerical models, such as the detailed radiative 
hydrodynamic models used here, together with detailed 
coronal models similar to those that have recently been 
successfully applied to the Sun (Gudiksen & Nordlund 
2001). By constructing a model that connects the surface 
and subsurface magnetic field to the corona, it is our hope 
that the energetic events observed in the coronae of M- 
dwarfs will provide a link to understanding the dynamo. 
This goal will be achieved by employing a sequence from a 
radiative hydrodynamic model as a lower velocity bound- 
ary condition for a coronal model, yielding an output that 
can be compared directly to observations. The unknown 
magnetic surface topology enters as a "free parameter" 
as an initial lower boundary condition that is driven to 
reconnect by the braiding surface motions. 
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